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Abstract. Using Monte Carlo simulations and published radial velocity surveys we have constrained the frequency and 
separation (a) distribution of very low mass star (VLM) and brown dwarf (BD) binary systems. We find that simple Gaussian 
extensions of the observed wide binary distribution, with a peak at 4au and 0.6 < o-iog{a/au) < 1-0, correctly reproduce 
the observed number of close binary systems, implying a close (a < 2.6 au) binary frequency of 17-30 per cent and overall 
frequency of 32-45 per cent. N-body models of the dynamical decay of unstable protostellar multiple systems are excluded 
with high confidence because they do not produce enough close binary VLMs/BDs. The large number of close binaries and 
high overall binary frequency are also completely inconsistent with published smoothed particle hydrodynamical modelling 
and argue against a dynamical origin for VLMs/BDs. 
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1. Introduction 

Very low-mass stars (VLMs, M < 0.15 M©) and brown 
dwarfs (BDs, M < 0.075 A/q) are now known to be common 
in our galaxy, possibly even outnumbering stars of higher 
mass. Yet their origins are still mysterious because the typical 
Jeans mass in a giant molecular cloud is an order of magni- 
tude more massive. 

Recent ideas, suggest that "turbulent" fragmentation 
might allow enhancements in local densities that can result 
in the formation of very small cores (e.g. Padoan & Nord- 
lund 2004). However, the biggest problem for the formation 
of VLMs and BDs may be how to prevent these very small 
cores subsequently accreting too much and becoming higher 
mass stars (Boss 2002). One solution is if VLMs and BDs 
were originally part of unstable multiple systems and, as the 
least massive components, are ejected on timescales too short 
for significant mass accretion (Boss 2001; Reipurth & Clarke 
2001). 

These ideas have proved very persuasive and it is possible 
that observations of accretion disks or VLM/BD kinematics 
and spatial distributions (see other contributions in this vol- 
ume) may ultimately provide crucial observational tests. A 
further observational constraint is offered by the frequency 
and separation distribution of binary systems and there have 
been many investigations along these lines in the last few 
years. 
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Binary systems should preserve information about their 
formation, particularly their dynamical origins. A plethora of 
studies have been published, which use high spatial resolu- 
tion imaging, either Hubble Space Telescope or adaptive op- 
tics in the near-infrared, to find the frequency of binary sys- 
tems among nearby field VLMS/BDs which have separations 
a > 1 au (e.g. Bouy et al. 2003; Burgasser et al. 2003; Close 
et al. 2003; Siegler et al. 2005). The lower cut-off to the sepa- 
ration range that can be probed is of course a result of limited 
angular resolution. The results of Close et al. are quoted be- 
low, but all of these studies are in reasonable agreement. The 
main findings are that the binary frequency for separations 
a > 2.6 au is 15 ± 7 per cent (see Fig. 1). This frequency is 
lower than found for G to M-type field stars (Duquennoy & 
Mayor 1991; Fischer & Marcy 1992), but these earlier stud- 
ies counted binaries over a much wider range of separations, 
including close, spectroscopic binary systems. 

Of more importance may be that all these studies failed 
to find binary VLMs/BDs with a > 15 au. This has been 
seized upon as supporting evidence for the "ejection" sce- 
nario where qualitatively, one expects fewer wide binary sys- 
tems with low binding energies to survive the ejection pro- 
cess. There are several problems with this: (i) some wider 
VLM/BD binary systems are starting to be found (e.g. Luh- 
man 2004; Phan-Bao et al. 2005); (ii) some theoretical mod- 
els that produce BDs by what could be termed "ejection" 
are also producing some BD-BD binaries with wide separa- 
tions (Bate & Bonnell 2005); (iii) there may yet be alternative 



3 RESULTS AND DISCUSSION 
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Fig. 1. The observed distribution of VLM/BD binaries as a 
function of log a (from Close et al. 2003). Note that the point 
at log(a/au) = is significantly affected by incompleteness. 
Also shown are the various model log a distributions used to 
simulate the VLM/BD binary population (see text). They are 
from Umbreit et al. (2005) or are Gaussian distributions with 
peaks at 4 au and with various indicated widths (truncated for 
a > 10 au). The model distributions are normalised such that 
the binary frequency is 15 per cent for 2.6 < a < 10 au. 



means of breaking up wide binaries during their first few Myr 
in a dense cluster environment and then one would not expect 
to see them in the field or older clusters. 

The frequency of close (a. < 2.6 au) binary systems may 
provide more secure evidence, but little progress has so far 
been made. Guenther & Wuchterl (2003) found 3 radial ve- 
locity (RV) variables among 24 field VLMs/BDs; Joergens 
(2005) found 2 RV variables from a sample of 1 1 VLMs/BDs 
in the Chamaeleon I star forming region; and Kenyon et al. 
(2005) found evidence for 4 close binaries in a sample of 
~ 60 VLMs/BDs in the young Sigma Ori cluster. Taken to- 
gether, these results hint at a surprisingly high (> 10 per cent) 
frequency for close VLM/BD binaries. In this contribution 
we examine this possibility with some statistical rigour and 
conoment on the imphcations of our results for the formation 
mechanism of VLMs/BDs. 



the value corresponds to a probability, p, of a constant RV 
that is less than 10^"^. This threshold was chosen to avoid 
contaminating what turns out to be a low number of binaries 
with any spurious statistical detections. During the course of 
this procedure we found, by investigating the distribution of 
p values, that the RV uncertainties were underestimated in 
the case of the Guenther & Wuchterl (2003) and Kenyon et 
al. (2005) data (for details see Maxted & Jeffries 2005). Er- 
rors amounting to 0.4 and 4.5 km s^^ respectively were added 
in quadrature to RV uncertainties from these two datasets. 
We identified four objects as binaries: 2MASSWJ21 13029- 
10094 and LHS292 from Guenther & Wuchterl, ChaHa 8 
from Joergens (2005), and KJNOM 72 from Kenyon et al. 

A Monte-Carlo simulation has been used to predict how 
many binary systems we ought to have found, given the sen- 
sitivity and time sampling of the data and for various assump- 
tions about the frequency, and distributions of a, eccentricity 
and mass ratio for binary systems. The details are provided 
by Maxted & Jeffries (2005) but briefly, we have assumed 
various Gaussian extensions of the observed log a distribu- 
tion of VLM/BD binary systems, peaking at 4 au and with 
widths ranging from (Tiog(a/au) = 0.6 to (Jiog(a/au) = 1-53 
(corresponding to the width of the field G-star distribution 
found by Duquermoy & Mayor 1991). These distributions 
were truncated at a = lOau in order to better match what 
is known about wider binary systems. We also trialled the 
prediction from a triple-body decay calculation for BD bina- 
ries from Umbreit et al. (2005) (see Fig. 1). The eccentricities 
were taken as e = for binary periods less than 10 days and 
then randomly distributed between < e < emax for longer 
periods, with emax = 0.6 or Cmax = 0.9. The mass ratio 
distribution was assumed either uniform for 0.2 < gr < 1.0 
(the "qflat" distribution) or uniform for 0.7 < < 1.0 (the 
"qpeak" distribution). These q distributions seem reasonable 
compared with what is already known about wider low-mass 
binary systems (e.g. see Bouy et al. 2003). It is possible, that, 
as in more massive binary systems, there is a shift towards 
equal masses among very close binary systems. Simulated 
targets were randomly assigned binary status on the basis of 
a total binary frequency which was normalised such that the 
frequency was 15 per cent for a > 2.6 au and hence matched 
the observational constraints for wider binaries. 



2. Analysis 

We have collected RV measurements for VLMs/BDs from 
three main sources in the literature. Consideration is re- 
stricted to those objects with RV measurements at > 2 
epochs. There are 24 field objects with 0.06 < M/Mq < 
0.10 from Guenflier & Wuchterl (2003), 10 objects wifli 
0.05 < M/Mq < 0.10 from Joergens (2005) and 14 objects 
with 0.045 < M/Mq < 0.11 from Kenyon et al. (2005). 
Masses were either estimated from spectral types (assuming 
an age of 1 Gyr for field stars) or from isochrones in colour 
magnitude diagrams. 

We identified possible binary systems by calculating the 
value of for a constant value fitted to the RV measure- 
ments of each object. Potential binaries are those for which 



3. Results and Discussion 

The outputs from our modelling are a prediction of the sen- 
sitivity of each of the target observations (for each star) to 
binaries of a given separation (averaged over all mass ra- 
tios and eccentricities). We plot the average sensitivity for 
all targets in Fig. 2 and also show how the data from each of 
the three datasets contributes to this average (i.e. the lower 
curves are the sum of the detection efficiencies for each tar- 
get in that datsaset, divided by the number of targets in all 
three datasets). It is clear from this that we are sensitive to 
binaries with a < 1 au. It is also clear that we cannot simply 
say binaries are all detectable out to some limiting separa- 
tion and cannot be detected beyond that. Even as a function 
of log a, the detection efficiency rolls over rather gradually, 



3 RES ULTS AND DISCUSSION 



1,0 






' ' 1 1 


' ' ' 1 


0.8 






¥ \\ 

\.Vl A ^^s^ 


qpeak, e<0,6 

qflat. e<0.6 

qflal, e<0.9 - 

JOS 

K05 

GW03 


0.6 










0,4 






V 




0,2 




--■v..,.--| 


■J \ ■ \ -"""^ 




0,0 






■ . 1 1 , , , — 





-3 -2-10 1 

log (o/ou) 

Fig. 2. The average binary detection efficiency. The upper 
three curves show the effects of differing models for the 
mass ratio and eccentricity distributions (see text). The lower 
curves show the contributions to the average from the each of 
the three RV surveys considered (i.e. the sum of these gives 
the average). To illustrate the binary population we might be 
sensitive to, we show the two extreme model a distributions 
we have investigated (see Fig. 1). 

and combined with the (probable) increasing number of bi- 
nary systems as log a increases, this demonstrates why a de- 
tailed simulation is demanded. From Fig. 2 we can see that 
narrower a distributions such as that of Umbreit et al. (2005) 
predict far fewer binaries will lie within our sensitive range. 
This is further illustrated in Fig. 3 where we show the simu- 
lated distribution of logp values (where p is the probability 
that a set of RV measurements represent a constant value) for 
a million binary systems with orbital separations taken from 
the various models. The number of detected binary systems 
(with logp < —3) is expected to vary by an order of magni- 
tude depending upon the adopted a distribution. 

Table 1. The number of binaries, nbin, predicted from the 
Monte-Carlo simulations. The total binary fraction is /b and 
probabilities (as percentages) are quoted for whether the sim- 
ulations would result in 4 or more, or 4 or fewer detected 
binaries in the observational dataset. 
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The key numerical results are given in Table 1, where we 
give the a distribution assumed, the implied total binary fre- 
quency, the range for the expected number of detected binary 
systems (summed over the three datasets for various combi- 
nations of the e and q distributions) and the probabilities that 
we would have observed either 4 or more, or 4 or fewer binary 
systems in the real data. We find that the chosen e or g distri- 
bution make little difference to the number of binary system 
detections expected (as implied by the similar average detec- 
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Fig. 3. The distribution of logp (averaged over all the objects 
included in the RV surveys and where p is the probability of 
the measured RVs being consistent with a constant value) for 
a simulated population of 10^ binaries. Separate curves are 
shown for binary separations which were drawn from each of 
the considered model a distributions. The straight line is the 
predicted locus for single stars. 

tion efficiency curves shown in Fig. 2). If binaries are really 
concentrated towards lower q values then we would expect 
to detect fewer of them and so would require a higher binary 
frequency to explain the data. Conversely, the minimum bi- 
nary frequency would arise from a population with q = 1. 
There is no evidence that either of these extremes are appro- 
priate. The upper end of the quoted range is produced by the 
"qpeak" q distribution. We have corrected the expected num- 
ber of binary systems for the well known bias in favour of 
binary detection in magnitude limited samples of field stars 
(for details see Maxted & Jeffries 2005). 

The expected number of detected binaries is very sensi- 
tive to the a distribution, ranging from < 1 for the Umbreit et 
al. (2005) distribution to 11 for the (Tiog(a/au) = 1-53 distri- 
bution. The narrowest (Umbreit et al.) a distribution is ruled 
out by the observations at the > 99 per cent level. There 
are too many binaries detected in the RV data. The broadest 
(o'iog(a/au) = 1-53) a distribution is also highly unlikely for 
the opposite reason. The intermediate width distributions are 
in satisfactory agreement with the observed number of close 
binaries and imply overall binary frequencies of 32^5 per 
cent, with 17-30 per cent at a < 2.6 au. This high close bi- 
nary frequency is larger than for either G-dwarfs (14 per cent, 
Duquennoy & Mayor 1991) or M-dwarfs (~ 10 per cent, Fis- 
cher & Marcy 1992). However, this is merely a continuation 
of the trend which is already apparent in the binaries with 
2.6 < a < 10 au. 

To clarify some points that might puzzle the reader, we 
note that: (i) We do not know the correct functional form 
for the a distribution, we have merely assumed a plausible 
form. To minimize the implied binary frequency requires an 
a distribution that is concentrated towards small a but which 
"joins on" to the observed distribution at large separations. 
Of course one could hypothesize a discontinuous distribution 
- for example a secondary peak of binaries at a < 0.1 au. 
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Whatever distribution is assumed however, the binary fre- 
quency for a < 2.6 au cannot be less than the 4/48 bina- 
ries we have observed even if we were 100 per cent efficient 
at detecting binaries out to that separation (which were are 
not by some margin - see Fig. 2). (ii) The model a distri- 
butions are normalised to 15 per cent for a > 2.6 au, but 
this number is uncertain. If it were increased to 22 per cent 
then tests show that narrower distributions would be allowed 
(though the Umbreit et al. distribution would still be ruled 
out at high significance), there would be slightly fewer bi- 
nary systems at a < 2.6 au, but the overall binary frequency 
would be very similar. Conversely, a binary frequency of 8 
per cent for a > 2.6 au would require broader a distribu- 
tions to explain the number of observed RV variables, more 
close binary systems, but again the overall binary frequency 
would be almost the same, (iii) The 17-30 per cent figure for 
close binaries is not a formal error range - it is the range of 
frequencies implied if the Gaussian forms assumed for the a 
distribution are correct and the observed binary frequency is 
4/48. (iv) q = I binaries may be harder to detect than we have 
calculated. Our simulations assume the secondary star is "in- 
visible", but clearly an equal mass binary system will have 
an SB2 spectrum. For close systems this will be obvious, but 
more separated systems, where the orbital velocity amplitude 
is smaller than the spectral resolution (certainly the case at 
a ~ 1 au for the datasets considered here), would not imme- 
diately reveal themselves and might appear to have constant 
RV. This effect can only act to increase the binjiry frequency 
inferred from the observations. 

Our results pose considerable difficulties for current ideas 
of VLM/BD formation. When multiple systems form by frag- 
mentation, the minimum separation should be a ^5 10 au due 
to the opacity limit. Closer binaries subsequently require the 
operation of dynamical or hydrodynamical hardening mecha- 
nisms (Bate, Bonnell & Bromm 2002). N-body decay models 
of unstable multiple systems produce very few close binary 
systems (e.g. Sterzik & Durisen 2003; Umbreit et al. 2005). 
This is because dynamical interactions capable of harden- 
ing a VLM/BD binary usually result in the ejection of one 
of the components and its replacement by a more massive 
object. Hydrodynamical hardening process are of course not 
included in pure N-body models. The high close binary fre- 
quency we have deduced indicates that these mechanisms 
must be important and effective. 

Smoothed particle hydrodynamic (SPH) models also have 
problems. Bate et al. (2002) and Bate & Bonnell (2005) 
present simulations which produce only ~ 8 per cent bina- 
rity among BDs. Whilst the rarity of wide separation systems 
is approximately reproduced, the simulations also yield very 
few binaries with a < 5 au. Even though the SPH models 
are currently hmited to resolutions of 1 au, the lack of sys- 
tems produced with separations of a few au indicates that 
they should not be common even when computational resolu- 
tions improve. It appears to be a common feature of both N- 
body and SPH models that the overall frequency of VLM/BD 
binary systems formed through the decay of initially unsta- 
ble multiple systems is much lower than in higher mass stars 
(Delgado-Donate et al. 2004; Hubber & Whitworth 2005). 



4. Conclusions 

From the detection of 4 close VLM/BD binary systems in 
several published RV surveys, we have been able to place 
strong constraints on the binary frequency and separation dis- 
tribution of such systems. We find that there are too many 
detected close binaries to be consistent with the rather nar- 
row separation distributions predicted by N-body models of 
dynamically decaying multiple systems. Instead we find that 
broader Gaussian distributions with 17-30 per cent binarity 
for a < 2.6 au are favoured and overall binary frequencies 
of 32^5 per cent. The total binary frequency for VLMs/BDs 
is therefore only marginally lower than (if not similar to) that 
for M-dwarfs, but significantly lower than for G-dwarfs. Con- 
sidering only separations a < 10 au, the binary frequency of 
VLMs/BDs is higher than for either G- or M-dwarfs. 

The neglect of hydrodynamical hardening processes in 
N-body models may be responsible for their lack of pre- 
dicted close binaries. However, recent SPH models also do 
not predict a significant population of very close BD binaries 
and yield overall binary frequencies of only ~ 8 per cent. 
The high overall binary frequency and the frequency of ob- 
served close binary VLMs/BDs do not favour the "ejection" 
hypothesis or similar models that produce BDs from the de- 
cay of dynamically unstable multiple systems. A mechanism 
is required that hardens VLM/BD binaries without disrupting 
them. 

At present, these conclusions are hmited only by the 
small numbers of detected binaries in the RV datasets. It is 
now important to increase these samples with further high 
precision RV surveys of VLMs/BDs. Such surveys will re- 
quire sparse monitoring of large samples with baselines of 
a year or more and RV precisions of ~ 1 kms^^. Not only 
can the techniques we have demonstrated yield the binary 
frequency, but modelling the observed logp distribution (as 
shown in Fig. 3) offers a means of determining the actual 
form of the a distribution at small separations, providing suf- 
ficient binaries can be detected. 
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